Theoretical study of doped-Tl2Mn2 07 and Tl2Mn2 07 under pressure 
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Using first-principles density functional based calculations, we study the effect of doping and pres- 
sure in manganese based pyrochlore compound, Tl2Mn207 that exhibits colossal magneto-resistive 
behavior. The theoretical study is motivated by the counter-intuitive experimental observation of 
suppression of ferromagnetic transition temperature upon apphcation of pressure and its enhance- 
ment upon substitution of Mn by moderate amount of nonmagnetic Sb ion. We also attempt to 
resolve the issue related to crystal structure changes that may occur upon application of pressure. 



I. INTRODUCTION 



Manganese-based pyrochlore compounds such as 
Tl2Mn207, form an interesting class of compounds which 
despite having similar colossal magneto-resistive (CMR) 
effect exhibit markedly different features compared to 
that of manganites. In contrast to manganites, these 
compounds are not mixed- valent, do not have any appre- 
ciable Jahn- Teller distortion^ and seemingly show metal- 
lic behavior'^ both below and above the magnetic transi- 
tion temperature Tc (« 140 K). Concerning the metallic 
behavior above T^ though, it is not totally clear whether 
it is a good metal or more like a bad metal or an insulator- 
like compound which arises due to difficulties related to 
sample preparation and defects present in the paramag- 
netic phase. Though the original data by Shimakawa 
et. al.^ suggests the metallic behavior above Tc, the 
optical conductivity data- and and photoemission and 
x-ray absorption^ near edge spectroscopy point towards 
an insulator-like behavior in the paramagnetic regime. 
These non-manganite like features rule out the gener- 
ally applicable double exchange mechanism or polaronic 
mechanism as the driving mechanism of ferromagnetism 
in this class of compound, and remained a mystery for 
quite some time. In absence of double exchange mech- 
anism, a strong contestant for the driving mechanism is 
the Mn-O-Mn mediated superexchange mechanism, par- 
ticularly since Mn-O-Mn form an angle^ of 133°, signif- 
icantly smaller than 180°. This is however contradicted 
by the experimental findings that the ferromagnetic Tc 
gets suppressed by application of pressure^ and gets en- 
hanced by moderate amount of introduction of nonmag- 
netic ion like Sb in place of Mn^. Employing N-th or- 
der muffin tin orbital (MTO) based NMTO-downfolding 
technique^, we have recently established^ that driving 
mechanism of ferromagnetism in this class of compounds 
is neither super-exchange or double-exchange, nor any 
exotic mechanism with nearest neighbor (NN) antifer- 
romagnetic (AFM) coupling dominated by long-ranged 
ferromagnetic (FM) interaction^^, but a kinetic energy 
assisted mechanism similar to that of double perovskite 
compounds like Sr2FeMo06, Sr2FeReOgii'i^ and Mn- 
doped GaAsi^. 



In the present paper, we give a detail account of the 
doping effect and the pressure effect from the electronic 
structure point of view, which lend further support to the 
existence of the proposed kinetic energy driven mech- 
anism in these compounds. While the changes that 
happen in the crystal structure upon doping is clear, 
the situation is rather ambiguous in case of application 
of pressure. The initial measurement hinted to a de- 
crease in Mn-O-Mn bond angle^ii^ alongwith the com- 
pression of Mn-0 bondlength. However the subsequent 
measurementri5 indicated an increase in Mn-O-Mn bond 
angle which was used in terms of super exchange mech- 
anism to explain the decrease in Tc. In view of these 
conflicting claims, we also take the opportunity to theo- 
retically address the effect of pressure on the structure of 
Tl2Mn207 and its implication on magnetism by detailed 
electronic structure calculations together with structural 
relaxations under pressure. 

The rest of the paper is organized in the following way. 
Section II deals with various methodologies that have 
been employed to compute the crystallographic struc- 
ture (in case of Tl2Mn2 07 under pressure) and to com- 
pute and analyze the subsequent electronic structure of 
Tl2Mn207 under doping and pressure. In section III, 
we present our results. This section is divided into sev- 
eral subsections. In section IIIA, the electronic structure 
of pristine Tl2Mn207 is presented and the proposed ki- 
netic energy driven mechanism is introduced as a refer- 
ence, which will form the basis of discussion in the follow- 
ing. Section IIIB is devoted to Tl2Mn207 under pressure 
, while section IIIC is devoted to Sb-doped Tl2Mn207. Fi- 
nally, we conclude the paper with a summary presented 
in section IV. 



II. METHODOLOGY 

We have investigated the electronic structure of 
pure Tl2Mn207, Tl2Mn2_:rSb:,07 and Tl2Mn207 under 
pressure by the tight-binding linear muffin-tin orbital 
(LMT0)>1& method, computed within the framework of 
the local spin density approximation (LSDA) of the den- 
sity functional theory (DFT). The employed basis set, 
consisted of Tl-6s, 6p and 5d states, Mn-4s, 4p and 3d 
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states, and 0-2s and 2p states. Three different classes of 
empty spheres were used to fill up the space. The self- 
consistent calculations were performed on a BZ mesh of 
12x12x12. 

The computed band structures were analyzed and in- 
terpreted in terms of NMTO-downfolding technique^. 
The NMTO technique, introduced and implemented in 
recent years, goes beyond the scope of standard LMTO 
technique by defining an energetically accurate basis with 
consistent description throughout the space involving 
both the MT spheres and the interstitial. As important 
feature of this technique is the energy-selective down- 
folding procedure where one integrates out degrees of 
freedom and keeps only few degrees of freedom active 
to define a few-orbital, low-energy Hamiltonian. The ef- 
fective NMTOs defining such Hamiltonian, serve as the 
Wannier or Wannier-like function corresponding to the 
selected low-energy bands. 

For the calculation of the exchange interaction 
strengths, J, a supercell involving a sixteen Mn atoms 
was constructed which amounts to enlarging the original 
unit cell by four times. This was achieved by replacing 
the original face centered cubic (FCC) lattice by the sim- 
ple cubic lattice with four times as many atoms as that 
in FCC. The basic structural unit of Mn sublattice is the 
Mn4 tetrahedra which share corners to form an infinite 
3D lattice. We calculated the LSDA total energies by 
considering different magnetic configurations of this sub- 
lattice. In total, six spin configurations have been chosen. 
One of the six configurations was FM and the other five 
were AFM. The details of the spin configurations chosen 
are as given in Table I in re^. The total energies are 
then mapped on to an effective Heisenberg Hamiltonian 
formed by the Mn^"* spins. The Hamiltonian considered 
till the 3NN(third nearest neighbor) interactions is: 



H 



Ji S^■SJ + J2 ^ S,.Sj + J3 ^ S,.Sj (1) 
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where Si denoted the spin-3/2 operator corresponding to 
the Mn spins at the site i, and Ji, J2 and J3 denote the 
NN, 2NN(second nearest neighbor) and 3NN magnetic 
exchange interaction strengths. Since we have restricted 
ourselves to only five AFM configurations amongst many 
possible AFM configurations, we computed seven differ- 
ent estimates of these J's along with their standard devia- 
tions employing set of three energy differences chosen out 
of total five different energy differences. This limitation 
can be overcome by taking as many AFM spin configu- 
rations as practically possible. The mean-field estimate 
of Tc is given by 



.^.f _ S{S + 1) 
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where Jq, is the net effective magnetic exchange interac- 
tion, given by 



and S — ^. ks is the Boltzman constant and 01,2:2, 23 
are the number of NN, 2NN and 3NN Mn pairs. For un- 
doped Tl2Mn207, zi — 6, Z2 — 12 and Z3 — 12. However, 
these numbers change in case of doped calculations where 
some of the Mn atoms gets replaced by the nonmagnetic 
Sb ion and therefore do not contribute in the calculation 
of magnetic energy. The application of mean-field the- 
ory is naturally expected to overestimate the T^,. Use of 
Monte Carlo simulation^- or Green Function method- 
can give better estimates and agreements with experi- 
mentally measured values. One needs to keep this in 
mind when comparisons are being made with the mea- 
sured Tc in the following sections. 

Tl2Mn207 occurs in cubic, face-centered lattice with 
Fd — 3m space group and two formula unit in the cell.— 
There are two types of oxygen, O' and O. While O oxy- 
gen provide the octahedral surrounding of the Mn atom 
which corner share to give rise to a geometry of a closed 
ring or cage, O' atoms form nearest neighbor bonds with 
Tl giving rise to TI2O' complex running through the cage 
like geometry as shown in left panel of Fig. 1. The 
structural parameters for the pristine Tl2Mn207, and Sb- 
doped Tl2Mn2 07 have been taken as the experimentally 
determined structure published in literature^^^^. How- 
ever, as mentioned already, the structural information 
for Tl2Mn2 07 under pressure is somewhat controversial 
which calls for structural optimization. For this purpose, 
we have adopted the pseudopotential method^° as imple- 
mented in the VASP—^ code. We have employed the pseu- 
dopotentials generated by the projector-augmented wave 
(PAW}2S, method. Wavefunctions have been expanded in 
terms of plane waves with a cut off energy of 400 eV. 
A A;-mesh of 6x6x6 was employed to perform the self- 
consistent calculations. In Tl2Mn2 07, Tl, Mn, O' and 
O atoms occupy the Wyckoff positions, 16d, 16c, 8b and 
48f respectively. Among these only the positions given by 
48f contain the internal parameter,?/. This parameter u 
defines the ZMn-O-Mn (see right panel of FigdJ. As u in- 
creases, the ZMn-O-Mn decreases and vice versa. When 
pressure is applied, it needs to be investigated whether 
the u parameter also varies along with the reduction in 
the bondlengths. If at all it changes, it has to be ascer- 
tained in what way it does so. In our calculation, the 
O position is relaxed by varying u, in order to find out 
the variation of total energy with angle formed between 
Mn, O and Mn atoms. In order to obtain the theoretical 
ground state, cohesive properties and estimate the pres- 
sure, Murnaghan22. equation of state has been employed 
to fit the LSDA total energies. The Murnaghan equation 
of state is given by 
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Jo = zi Ji -I- Z2J2 + Z3J3 



(3) 



where Vq is the equilibrium volume, Bq is the bulk mod- 
ulus and is given by Bq = —V{6P/6V)t evaluated at 
volume Vo- Bq is the pressure derivative of Bq also eval- 
uated at volume Vq. Bq provides a measure of stiffness 
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FIG. 1: (color on-line) Left panel: The Tl2Mn207 structure consisting of the ring formed by MnOe octahedra with Tl-0' rods 
passing through the holes of the rings. Right panel: Backbone of the structure realized by the ring like structure formed by the 
Mn-O-Mn network (the Tl-O' rod network is not shown for clarity). Gold, green and black or cyan spheres denote the Tl, Mn 
and O atoms. Black and cyan spheres correspond to oxygen (O) atoms with two different oxygen (O) position parameter u. 



of the material upon increasing pressure. Typical val- 
ues for the parameter Bq are between 4 and 7. Pressure, 
P, corresponding to volume, V, occupied by a cell of a 
particular lattice parameter is given by 



P{V) 



Bo 
B' 



[Vo/vr- - 1 



(5) 



III. RESULTS 

A. FM Electronic Structure of Pristine Tl2Mn207 
and the Underlying Mechanism 

The LSDA electronic structure of ferromagnetically or- 
dered TlaMnaOy, computed with TB-LMTO method is 
shown in Fig[2l This figure shows both the band struc- 
ture and the orbital projected density of states (DOS). 
Consistent with the nominal valence of and Mn'*+, 
oxygen dominated states are totally occupied, while the 
crystal field split Mn t2g states are occupied in the major- 
ity spin channel and close to empty in the minority spin 
channel. Mn eg states consistent with ?>d? configuration 
of Mn, remain empty in both the spin channels. Tl bands, 
consistent with nominal valence also remain almost 
empty. However, interestingly in the majority spin chan- 
nel, the Tl-6s state mixed with 0'-2p states overlap with 
Mn Cg manifold spanning an energy range of about 2-4 
eV— , while in the minority spin channel, they overlap 
with nearly empty Mn t2g dominated states, which give 
rise to highly dispersive band in bottom of this manifold. 



crossing the Fermi level, Ei?, as seen in lower, left panel 
of FigO This hints to a large spin-splitting at the Tl 
site. This is unusual in the sense, one would normally 
expect Tl to be essentially nonmagnetic with tiny spin 
splitting. This unusual feature plays the key role in pro- 
viding the understanding of the driving mechanism of 
ferromagnetism in Tl2Mn2 07. This has been explained 
in detail in ref®. As has been shown in re^, the unhy- 
bridized effective Tl-O' level is positioned in between the 
spin-split Mn t2g states and on turning on the hybridiza- 
tion, it induces a renormalized spin-splitting within the 
Tl-O' level which is directed in a opposite way compared 
to the spin-splitting at the Mn site. This happens due to 
pushing up of the Tl-O' spin-up state and pushing down 
of the Tl-O' spin-down states, driven by coupling with 
the Mn t2g states of the same symmetry. The energy gain 
contributed by this hybridization induced negative spin 
polarization of the otherwise nonmagnetic element, is the 
central concept in this novel mechanism. This mecha- 
nism ensures a specific spin orientation between the lo- 
calized Mn t2g spins and the mobile carrier formed by the 
hybridized renormalized state, which in turn aligns the 
spins in Mn sublattice. This is a general mechanism and 
found to be operative in case of a number of compounds 
like SraFeMoOg, SrsFeReOgiiii^, Mn-doped GaAs^^ and 
CuCraS^S^ etc. 

In the following subsections, we will present a detail ac- 
count of the changes in the electronic structure and the 
variation in Tc of Tl2Mn207 with doping and under pres- 
sure. We will argue that the counter-intuitive variation 
in Tc has its origin in the modified electronic structure 
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FIG. 2: (color on-line) Left panels show the spin-polarized 
band structure computed within LSDA for Tl2Mn207 . The 
zero of the energy is set at the LSDA Fermi energy, E f • Right 
panels give the spin-polarized total (black), Tl-6s (red), Mn- 
3d (green), 0-2p (orange) and 0'-2p (blue) density of states. 
Upper and the lower figures correspond to the majority and 
minority spin contributions respectively in both the panels. 
Bars on the right signify the energy range for which various 
characters dominate. 



and follows naturally from the framework of proposed 
hybridization induced mechanism of ferromagnetism as 
was established in ref^ and reviewed in the above. 



B. Tl2Mn2 07 under pressure 

As explained before, to investigate Tl2Mn207 under 
pressure one needs to first settle the issues related to 
crystal structure. To be precise, does applied pressure in- 
crease or decrease or not alter the ZMn-O-Mn ? For this 



purpose, we have carried out a series of structural op- 
timization calculations with pseudopotential mcthod^». 
At first the cohesive properties were calculated within 
LSDA. Following the recent claim^ that application of 
on-site U can also infiuence the cohesive properties of 
transition metal oxides, we have repeated our calcula- 
tions also with LSDA+U functional. 
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FIG. 3: Total energy vs. volume of the unitcell calculations 
within LSDA and LSDA-I-U. The symbols denote the DFT 
total energies and the solid and dashed lines depict the total 



energies obtained from fitting with Murnaghan 
state (EOS). 



equation of 



LSDA total energy vs volume calculation, shown in 
FigEl has a minimum corresponding to lattice constant 
of 9.881 A. Considering the experimental value of the lat- 
tice constant is 9.892 A, this gives a deviation of about 
0.1 %, which is within the limit of LSDA overbinding. 
The value of the bulk modulus, obtained by fitting the 
total energy vs volume curve with Murnagahan equa- 
tion of state^^ is 230.6 GPa, which is of the same order 
as the experimentally measured value of «200 GPa in 
Tli.8Cdo.2Mn207i^. Bq is -1-4.84, much higher than that 
in Cd doped Tl2Mn207 but within the expected limit. 
Applying an on-site Coulomb interaction of U=2 eV and 
J=0.7 eV on the Mn site leads to a slightly better agree- 
ment with the experimental values. The lattice param- 
eter changes to 9.888 A, Bq and Bq remain unaltered. 
In view of the slightly improved results of LSDA-I-U, we 
have used the Bq and Bg as estimated by LSDA+U cal- 
culations to compute the pressure as given by Eq.(5). 

Till now we have considered the experimental ZMn-O- 
Mn as the one for the theoretical ground state. One needs 
to relax this angle in order to obtain the true theoretically 
predicted ground state. The co-ordinate of O atom in the 
basis is given by (u,. 125,. 125). The u parameter will push 
O atom in or out depending on its value, thereby chang- 
ing the angle formed between Mn, O and Mn atoms. The 
u parameter is relaxed in the subsequent calculations, 
and the LSDA (for 9.881 A) and LSDA-hU total energies 
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T T\ A 


LSDA+U EXPT. 


a(A) 


9.881 


9.888 


9.892 


u 


0.3214 


0.3227 


0.3254 


Mn-0 (A) 


1.886 


1.890 


1.901 


ZMn - O - Mn(deg 


) 136.01 


135.29 


133.8 


Bo (GPa) 


230.6 


230.7 




b; 


+4.84 


+4.85 





154 



TABLE I: Calculated ground state structural parameters of 
Tl2Mn207 within LSDA and LSDA+U compared with the 
experimenta]^ data at ambient condition. 



(for 9.888 A) are computed varying the u parameter. The 
minimum in the E vs m curve (see Fig|3]) is obtained for 
u=0.3214 and 0.3227 respectively. The ZMn-O-Mn in the 
ground state obtained within LSDA and LSDA+U frame- 
work are given in Table I along with the corresponding 
experimental values. Slightly better results are obtained 
for the LSDA+U calculations, a trend already observed 
in calculation for equilibrium lattice constant. The result 
match with experimental values upto the second decimal 
point for u and within 1.5° for Mn-O-Mn angle. 
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FIG. 4: Relaxation of the Oxygen-0 position parameter u at 
ambient condition within LSDA and LSDA+U in Tl2Mn2 07. 
The symbols are the DFT total energies and the solid and 
dashed lines are cubic spline interpolated. 

Having established the ground state, in the next step 
we relax the O position for lattice parameters a=9.846 
and 9.69 A (LSDA+U) which amount to applying a pres- 
sure of 3.04 GPa and 16.28 GPa. This wiU clarify whether 
the Mn-O-Mn angle also changes with the decrease in the 
interatomic distances. The minimum obtained from the 
total energy vs. u plots, a shown in Fig. [Slare listed in Ta- 
ble IL We find that on application of pressure, the ZMn- 
O-Mn reduces insignificantly, by about 0.1° and 0.4° 
for pressure of 3.04 GPa and 16.28 GPa respectively. A 
reduction of the angle by about 0.6° occurs for the same 
applied pressure (16.28 GPa) within LSDA. We thereby 
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FIG. 5: Total energy vs u, i. e. relaxation of the Oxygen (O) 
position parameter u at ambient condition and under pressure 
of 3.04 and 16.28 GPa within LSDA+U in Tl2Mn207. The 
symbols are the DFT total energies and the solid, dashed and 
dashed-dotted lines are cubic spline interpolated. 





AMBIENT 3.04 GPa 


16.28 GPa 


u 


0.3227 0.3230 


0.3234 


Mn-O(A) 


1.890 1.883 


1.854 


IMn - O - Mn (deg 


) 135.29 135.18 


134.90 



TABLE II: Structural parameters obtained within LSDA+U 
at ambient and two different pressures upon relaxation of the 
Oxygen(O) position parameter u in Tl2Mn207. 



establish that, regardless of kind of DFT one is using, on 
applying pressure ZMn-O-Mn does not increase rather 
decreases though almost insignificantly. The major ef- 
fects comes from reduction in the bondlengths. In the 
following calculations, we have therefore kept the ZMn- 
O-Mn unchanged upon application of pressure. 

We have carried out LSDA total energy calculations for 
various spin arrangements of Mn atoms as explained in 
section II. The relative LSDA total energies, as given in 
Table III are mapped on to Eq. (1) in order to extract the 
J's. The magnetic exchange interaction strengths are cal- 
culated and shown in Table IV for almost 2 % reduction 
in experimental lattice parameter, amounting to a pres- 
sure of 16.28 GPa. Considering pressure of 16.28 GPa, 
we observe a substantial decrease in the NN FM inter- 
action. The 2 and 3 NN interactions remain almost the 
same. This leads to a mean-field Tc of about 140 K, lead- 
ing to a decrease of about 41 K. The gradient is about 
-2.5 K/GPa. The gradient calculated from another set 
of calculations at 3.04 GPa is about -2.3 K/GPa. One 
can approximate a linear relation between T^ and ap- 
plied pressure. Thus the average gradient is about -2.4 
K/GPa. This is not an unreasonable estimate consider- 
ing the mean-field overestimation. The ratio of exact Tc 
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AE(meV) 




Normal 


3.04 GPa 


16.28 GPa 


FM 











AFMl 


9.55 


9.18 


7.36 


AFM2 


13.98 


14.17 


13.05 


AFM3 


12.24 


11.49 


8.82 


AFM4 


17.59 


22.5 


16.24 


AFM5 


18.39 


17.7 


14.04 



TABLE III: Relative LSDA energies per Mn ion in meV in FM 
and five AFM spin configurations for Tl2Mn207 at normal, 
3.04 GPa and 16.28 GPa pressures. All energies are converged 
upto 0.01 meV per Mn ion. 

and T^'f is estimated as 0.79 and 0.81 for BCC and FCC 
lattice with coordination 8 and 12 respectively. Consid- 
ering the average coordination of 10 for Mn sublattice 
in Tl2Mn207, one would expect the gradient in Tc to 
be about -1.9 K/GPa. This is a bit larger than the ex- 
perimental values, which are -1.6 K/GPcfc^ and about -1 
K/GPa^^ . A plausible explanation for this somewhat 
overestimation may be due to computational limitations 
imposed by the number of spin configurations. A better 
estimate of the gradient might have been reached had 
all the possible spin configurations been considered for 
calculation of the J's. Therefore one can conclude that 
within the limitation of our method, we have established 
the trend in ferromagnetism of Tl2Mn207 with applied 
pressure correctly. 



Pressure Appl. 


Ji 


h 


Js 


Jo 


rpm/ 

^ c 


(GPa) 


(meV) 


(meV) 


(meV) 


(meV) 


(K) 





-2.52 


-0.11 


-1-0.33 


12.48 


181 


3.04 


-2.48 


-0.16 


4-0.40 


12.06 


174 


16.28 


-1.99 


-0.16 


-F0.35 


9.66 


140 



TABLE IV: NN, 2NN, 3NN and effective exchange interac- 
tion strengths and mean-field Curie Temperature calculated 
within LSDA at ambient condition and at 3.04 and 16.28 GPa 
for Tl2Mn2 07. 

In FiglB] we show the density of states of Tl2Mn207 
under 16.28 GPa pressure alongwith that of Tl2Mn207 
under ambient pressure for comparison. Upon applica- 
tion of pressure, the bondlengths shorten which enhances 
the Mn-0 and Tl-0' interaction and as a result net band- 
width expands and t2g-eg splitting enhances. However 
the key effect is seen in the inset, where we show the 
Mn-d and Tl-s projected density of states in an energy 
window close to E^? in solid and dashed lines respectively. 
They are shown in two different colors, corresponding to 
ambient and 16.28 GPa pressure. We notice that the 
relative proportion of Tl-0' in the hybridized, dispersive 
band has reduced significantly upon application of pres- 
sure. It has reduced from 45 % in case of ambient pres- 
sure to 36 % in case of 16.28 GPa pressure. This indicates 
reduction of Mn-Tl-0' hybridization upon application of 
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FIG. 6: (color on-line) Spin-polarized density of states at am- 
bient condition (black) and 16.28 GPa (red) calculated within 
LSDA for Tl2Mn207. The upper (lower) panel corresponds 
to the majority (minority) spin. The negative of density of 
states has been plotted in the lower panel for clarity. Inset 
shows minority spin channel density of states for the energy 
region near Fermi level at ambient and 16.28 GPa. The solid 
and dashed lines show the contribution from Mn - 3d and 
Tl-6s-|-0'-2p states respectively within LSDA. 



pressure. Since the underlying mechanism of ferromag- 
netism is based on the Mn-Tl-0' hybridization, weaken- 
ing of such hybridization would cause a decrease in T^ as 
supported by the detail calculation presented above. The 
change in the hybridization is proportional to the change 
in the hopping interaction strength, t and inversely to 
the change in the energy level separation A. Applica- 
tion of pressure causes shortening of bondlength, hence 
one would expect hopping to increase and in turn the 
hybridization to increase, opposite to what is observed 
in the calculation. For this we performed the NMTO- 
downfolding calculation keeping Mn-3d, Tl-6s and 0-2p 
degrees of freedom active and integrating out the rest, 
including O' degrees of freedom to define a Tl-0' effec- 
tive level. The real-space representation of that down- 
folded Hamiltonian in the NMTO-Wannier function ba- 
sis for Tl2Mn207 under application of 16.28 GPa pres- 
sure shows the Tl-0' effective level to shift down by 0.3 
eV compared to that under ambient pressure. This in- 
creases the A and substantially reduces the Tl-O'-Mn 
hybridization, substantial enough to compensate the in- 
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crease in tsda- hopping interaction between Mn-t2g and 
the effective Tl-s orbital and to produce a net negative 
contribution to the hybridization. This is an interesting, 
counter-intuitive situation which can only be unraveled 
by a detail analysis of the underlying electronic structure 
as has been done in the present study. 



C. Sb-doped Tl2Mn207 

To mimic the doping effect, we have carried out cal- 
culations using supercell containing sixteen Mn atoms. 
Construction of such supercell has been described in sec- 
tion II. Out of these, sixteen Mn atoms, some are sub- 
stituted by Sb to simulate the effect of doping. There 
are two doping levels that have been investigated. In one 
case, one of the Mn atom is substituted by an Sb atom, 
and in another case, two of the Mn atoms are replaced 
by two Sb atoms. The former amounts to doping level 
of x=0.125 or 6.25% and the latter amounts to doping of 
x=0.25 or 12.5 %. These are not exactly the doping lev- 
els that has been investigated experimentally but close 
to them {i.e. x = 0.1 and 0.2 respectively). 

Substitution of Mn''+ by larger cation Sb''^+ does not 
cause any significant changes in the crystal structure, 
other than the fact that lattice expands slightly keep- 
ing the various angles unchanged. The lattice parame- 
ters corresponding to doping levels of x=0.125 and 0.25 
are obtained by interpolation of the experimentally de- 
termined lattice constants at various doping as shown in 
the inset of Fig.l. in refi. The lattice parameters for x 
=0.125 and 0.25 are increased by 0.3% and 0.9% respec- 
tively, compared to undoped lattice constant. 

To estimate the exchange interaction strengths for 
Tl2Mn2-2:Sb:r07 with x=0.125 and 0.25 we have carried 
out LSDA total energy calculations for various spin ar- 
rangements of Mn atoms as explained in section II. The 
relative LSDA total energies, as given in Table V are 
mapped on to Eq. (1) in order to extract the J's. At 6.25 





AE (meV) 




x=0 


x=0.125 


x=0.25 


FM 











AFMl 


9.55 


16.50 


23.96 


AFM2 


13.98 


26.07 


29.39 


AFM3 


12.24 


23.18 


31.53 


AFM4 


17.59 


33.30 


35.69 


AFM5 


18.39 


34.39 


40.11 



TABLE V: Relative LSDA energies per Mn ion in meV in FM 
and five AFM spin configurations for doped Tl2Mn207. All 
energies are converged upto 0.01 meV per Mn ion. 

% doping the average coordination numbers change to 
5.6, 11.2 and 11.2 respectively. In order to check whether 
the chosen disordered configuration has any influence on 
the trend, we have repeated the calculation with two Sb 



atoms in the unit cell for two cases. In one case, two Sb 
atoms were chosen as NN pairs while in the second case 
they have been chosen as 3NN pairs. The respective co- 
ordination numbers at 12.5 % doping are 5.14, 10.28 and 
10.86 for 3NN substitution and 5.28, 10.28 and 10.28 for 
NN substitution. In all the three cases, i.e for x=0.125, 
0.25 (NN) and 0.25 (3NN), the ground state is found to 
be FM. 



doping 


Ji 


J2 J3 Jo 


D(Ef) 


rpmf 
^ c 


cone. 


(meV) ( 


meV) (meV) (meV) (States/eV) 


(K) 


x=0.0 


-2.52 


-0.11 -hO.33 -12.48 


0.2 


181 


x=0.125 


-3.65 


-0.12 -0.41 -26.38 


1.0 


382 


x=0.25(3NN) 


-2.72 


-0.57 -0.85 -29.08 


1.4 


422 


x=0.25(NN) 


-2.74 


-0.56 -0.87 -29.30 


1.44 


423 



TABLE VI: Exchange interaction strengths Ji (NN), J2 
(2NN) and J3 (3NN) and effective exchange coupling strength 
Jo in meV, total density of states at Fermi level, D(Ef) 
and mean- field Curie Temperature, T^f for Tl2Mn2-^Sb^07, 
x=0.0, 0.125 and 0.25. 



Table VI gives the magnetic exchange interaction 
strengths for all the three cases and their corresponding 
Tc's and density of states at Ep, D(Ei?). For compar- 
ison, we also quote the numbers corresponding to pure 
Tl2Mn207 taken from ref^. Pristine Tl2Mn207 has a 
very strong NN FM interaction, a weak 2NN FM inter- 
action and a little stronger 3NN AFM interaction. With 
the introduction of doping, all the interactions became 
FM and their strengths also increase. On further in- 
crease of doping, the 2nd and 3NN FM interactions be- 
come stronger with a small decrease in the NN FM in- 
teraction. It is worth mentioning here that the results 
for x=0.25 doping with two NN Mn atoms substituted, 
does not change significantly compared to the result when 
the two Mn ions substituted were 3NNs, proving that 
choice of disordered configuration does not alter the gen- 
eral trend. Computation of Tc, using the values of Ji, 
J2 and J3, listed in Table I, show rapid increase in T^ 
on the initial doping which then attains a kind of satura- 
tion on further doping, considering upto moderate level 
of doping. This trend is in good agreement with experi- 
mental findings, but the quantitative values differ. This 
happens primarily because of the clustering effect not 
included in our calculation. To quote Alonso et. al^, 
" For low doping levels, we cannot disregard some type 
of magnetic clustering with almost pure Tl2Mn2 O7 (Tc — 
135K) and Sb rich regions". Such effects are bound to re- 
duce the global Tc measured experimentally. Also lattice 
relaxation effect could be important. Although experi- 
mentally measured crystal structure data shows the same 
space group symmetry as that of the undoped one with 
little variation in the internal u parameter, the substitu- 
tion of Mn ion by a larger cation will probably change 
the structure locally. Therefore structural relaxations of 
the doped compounds, which is beyond the scope of our 
present computational capacity, are needed to resolve the 
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issue completely. 
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FIG. 7: (color on-line) Spin-decomposed density of states 
for x=0.0 (black), 0.125 (red) and 0.25 (blue) doping in 
Tl2Mn2-a;Sbi07. Zero of the energy is set at the Fermi level 
of the undoped Tl2Mn207 compound. The upper (lower) 
panel corresponds to the majority (minority) spin. The neg- 
ative of density of states has been plotted in the lower panel 
for clarity. Inset shows the density of states in minority spin 
channel at the same doping levels for the energy region close 
to Fermi level. The dashed lines in the inset correspond to 
the respective Fermi levels. 

One needs to examine the electronic structure upon 
doping to understand the microscopic origin of the strik- 
ing effect described above. In Fig. [7] we show the density 
of states corresponding to Sb-doped Tl2Mn2 07 for dop- 
ing level of x=0.125 and 0.25. To appreciate the changes 
in the electronic structure upon doping, we also show 
the density of states for pure Tl2Mn207. Comparing the 
various density of states, we observe slight shrinking of 
the overall band-width due to the small expansion of the 
lattice. The t2g-eg sphtting at the Mn site decreases due 
to expansion of the Mn-0 bondlength. Smearing of the 
DOS, compared to pure case is also observed due to the 
effect of disordering. The Sb ions in the 5-1- state have 
completely filled 4d shell which are essentially core-like 
while 5s and 5p states remain more or less empty with 
no contribution to the states close to E^?. The most im- 
portant effect is that of electron doping. Replacement 
of Mn^+ by Sb^"*" introduces additional electrons in the 
system which have no other option but to populate the 



highly dispersive Tl-0'-Mn-t2g hybrid state in the mi- 
nority spin channel. Because of the highly dispersive 
nature of this band, the density of state in the minor- 
ity spin channel close to Ej? rises fast after the tail part 
which then attains a kind of small plateau before attain- 
ing highly peaked values corresponding to weakly disper- 
sive Mn t2g dominated states between 1-2 eV. This is 
shown in the inset of FiglT] With introduction of first 
Sb atom, the doped extra electrons while populating the 
dispersive minority spin band increases D{F,p) substan- 
tially. Upon increasing the number of Sb atoms from one 
to two, the Fermi level shifts towards right to accommo- 
date more electrons, but since it has almost crossed the 
highly dispersive region and falls within the plateau-like 
structure, D(Ep) does not change significantly. This is 
evident by quantitative estimates of D(Ei?) shown in Ta- 
bic I. In FiglXl we have shown only the case of x=0.25 
(NN), but the DOS corresponding to x=0.25 (3NN) is not 
very different and D{F,p) changes only marginally. Fol- 
lowing the perturbative treatment in terms of the transfer 
integraP'', t, between the magnetic and nonmagnetic site, 
the exchange coupling, J, within the proposed kinetic en- 
ergy driven scheme is given by, t^D{Ep) / A'^ , where A is 
the average energy separation between the magnetic and 
nonmagnetic level and D(Ei?) is the density of states at 
Fermi level introduced in the above. The real space rep- 
resentation of NMTO-downfolded Hamiltonian keeping 
Tl-6s, Mn -3d and 0-2p states active and downfolding 
the rest shows that the Tl-Mn hopping integral, t's and 
their energy level separation changes only slightly for 0.3 
% and 0.9% expansion of the lattice upon x=0.125 and 
0.25 doping. The major changes come from the sharp 
increase iiiiD(Eip) which explains the enhancement of J's 
and Tc's in turn. 



IV. CONCLUSIONS 



We presented a detailed study of Tl2Mn2 07 under 
pressure and Sb-doped Tl2Mn2 07 using first-principles 
electronic structure calculations. The analysis of the 
computed electronic structure shows that the counter- 
intuitive experimental observation of suppression of fer- 
romagnetic Tc upon application of pressure and its en- 
hancement upon moderate amount of doping by Sb in 
Mn sublatticc is naturally explained in terms of the hy- 
bridization induced mechanism proposed earlier—. Ap- 
plication of pressure reduces the hybridization between 
localized Mn t2g level and the delocalized Tl-0 effec- 
tive state due to the shifting of the position of the Tl-0 
effective level, which weakens the driving mechanism of 
ferromagnetism and thereby reduces the strength of ex- 
change coupling and the Tc. 
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In case of doping, on the other hand, the enhancement 
of Tc is induced by the charge-carrier doping of the sys- 
tem by Sb. This enhances the value of density of states 
at the Fermi energy significantly, thereby enhancing the 
exchange coupling and the Tc. Our conclusions were sub- 
stantiated by the explicit calculation of the J's and the 
Tc's. 

In case of Tl2Mn2 07 under pressure , we resolved the 
issue related to the change in Mn-O-Mn angle. We car- 
ried out structural relaxation by means of accurate pseu- 
dopotential based total energy calculations. Our theoret- 
ical results show, as opposed to the claim by Velasco et. 
ali^, the Mn-O-Mn angle decreases in agreement with 
the original suggestion by Sushko et. al.^ . However the 



decrease is only marginal and for all practical purposes, 
it may be assumed to remain same as that of Tl2Mn207 
in ambient condition. This rules out the use of angle 
variation argument to explain the reduction of Tc. 
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